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1.0  OVERVIEW  AND  SUMMARY  OP  USER'S  MANUAL 


1.1  INTRODUCTION  AND  GENERAL  COMMENTS 

;  This  report  was  prepared  as  an  introduction  to  using  the 
STAGSC-1  computer  program.  It  is  not  an  attempt  to  rewrite  the 
current  STAGSC-1  User's  Manual  [1];  however,  as  with  any  manual, 
there  are  certain  aspects  which  sometimes  are  not  easy  to  under¬ 
stand.  This  report  is  especially  addressed  to  the  user  that 
might  not  be  as  experienced  as  the  STAGSC-1  manual  presumes. 

Clarification  of  selected  input  data  cards  which  this 
reviewer  found  somewhat  confusing  are  covered  in  Section  1.2. 
Section  1.3  of  this  report  discusses  selected  output  messages, 
nomenclature,  and  options  which  will  aid  in  understanding  the 
STAGSC-1  output.  Section  2.0,  which  comprises  the  bulk  of  this 
report,  is  an  assembly  of  various  sample  runs.  The  input  data 
for  the  sample  problems  are  included  in  the  main  text  of  this 
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1.1.1  tfse  of  This  Report  J 

This  report  is  aimed  at  the  beginning  user  of  STAGSC-1  and 
possibly  the  beginner  at  structural  and/or  finite  element  analy¬ 
sis.  The  following  procedure  for  understanding  the  use  of  STAGSC-1 
is  recommended.  First,  review  this  report  without  reading  in 
detail.  Especially,  review  the  sample  problems  and  input  data 
format  presented  in  Section  2.  Next  read  Sections  1,  2  and  the 
first  few  pages  of  Section  3  from  the  STAGSC-1  User's  Manual  [1]. 

Next,  examine  the  first  example  (Section  2.1.1)  of  this 
report  and  look  up  any  card  in  Section  3  of  the  User's  Manual 
which  is  not  understood.  If  the  explanation  in  the  User's  Manual 
seems  confusing,  check  Section  1.3  of  this  report  to  see  if  addi¬ 
tional  clarification  of  the  card  has  been  covered.  Finally,  run 
the  sample  problem  and  examine  the  output.  Refer  to  Sections 
1.3.1,  1.3.2,  2.1.6,  and  Appendix  E  if  questions  concerning  the 
output  arise.  If  the  reader  Is  interested  in  the  complete  ouput 
of  the  ten  sample  problems,  he  should  contact  A1WAL/FIBRA. 


/  '  > 
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After  running  through  a  sample  problem,  one  can  proceed  in 
whatever  fashion  suits  one  best.  The  purpose  of  the  additional 
samples  in  Section  2  is  to  illustrate  various  capabilities  of 
STAQSC-1. 

It  should  be  mentioned  that  although  the  STAGSC-1  User’s 
Manual  does  not  contain  an  index  as  such,  the  mini-manual 
(Section  9)  seems  to  serve  the  same  purpose.  It  is  very  useful 
for  locating  the  proper  cards  without  rummaging  through  the 
voluminous  data  card  descriptions  (Section  3). 

1.1.2  Overview  of  the  STAGSC-1  Program 

The  STAGSC-1  program  originated  as  a  finite-difference  based 
structural  analysis  program.  Eventually,  it  was  converted  to  the 
present  configuration  which  is  entirely  finite  element  based. 
However,  some  of  the  nomenclature  which  was  used  in  the  generation 
of  meshes  and  grids  was  retained.  The  terms  "rows"  and  "columns" 
are  examples. 

The  STAGSC-1  computer  program  is  comprised  of  four  modules : 
STAGS1 ;  STAGS2;  P03TP;  and  STAPL.  The  STAGS1  module  is  a  pre¬ 
processing  module  which  performes  model  generation,  and  various 
preliminary  calculations.  Execution  of  this  module  typically 
precedes  the  execution  of  the  other  modules.  The  STAGS2  module 
performs  the  bulk  of  the  numerical  computations  associated  with 
finite  element  analysis.  STAGS2  module  performs  matrix  decompo¬ 
sition,  linear  and  nonlinear  stress  analyses,  eigenvalue  analyses, 
etc.  The  POSTP  module  is  a  postprocessing  module  for  determining 
secondary  solutions  from  previously  calculated  displacement  solu¬ 
tions  which  have  been  saved  on  the  restart  file  TAPE22.  The 
STAPL  module  is  a  module  for  plotting  undeformed  and  deformed 
geometries  and  solution  contour  plots.  This  manual  is  directed 
toward  the  input  and  output  data  associated  with  the  STAG31  and 
3TAC32  modules.  For  information  on  the  use  of  POSTP  and  STAPL, 
the  reader  is  referred  to  Section  5  of  [1], 
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STAGS-1  employs  unique  ’’substructures"  known  as  shell  units 
and  element  units.  A  shell  unit  is  a  collection  of  finite  ele¬ 
ments  and  nodes  which  may  be  defined  using  standard  shell  shapes 
or  with  user-written  subroutines.  The  user  defines  a  shell  unit 
by  specifying  limited  data  concerning  the  shell  and  letting 
STAGSC-1  generate  the  actual  internal  nodes  and  elements.  The 
element  unit  is  more  like  a  regular  finite  element  mesh.  Only 
one  element  unit  may  be  specified,  and  it  is  defined  by  node 
locations  and  element  connectivities  as  in  most  standard  finite 
element  codes.  The  one  element  unit  can  be  used  to  specify  any 
number  of  finite  elements  and  nodes,  their  locations  and  con¬ 
nectivities  to  other  shell  units  being  arbitrary.  In  the  user's 
manual  the  "A  through  R"  cards  are  concerned  mostly  with  defining 
shell  units,  while  the  "S  through  V"  cards  are  used  for  defining 
element  units.  There  are,  however,  certain  cards  in  the  "A 
through  R"  sequence  which  must  always  be  present,  regardless  of 
whether  shell  or  element  units  are  being  employed.  Paging  through 
the  user's  manual  will  allow  identification  of  these  cards. 

Another  apparently  unique  aspect  of  STAGSC-1  (although  it  is 
mostly  terminology)  is  the  use  of  "surface  coordinates."  These 
are  discussed  in  Section  2  and  are  referred  to  in  Figure  3.9  of 
[1].  Surface  coordinates  are  simply  alternative  coordinates  used 
to  describe  a  surface.  By  definition,  all  surfaces  can  be  des¬ 
cribed  by  two  independent  coordinates.  These  are  denoted  "surface 
coordinates"  in  STAGS.  For  example,  in  normal  nomenclature,  a 
cylinder  can  be  described  by  three  Cartesian  coordinates,  x,  y, 
and  z,  with  an  equation  relating  them,  4:hat  is: 

2  2  2  2 
x  +  y  +  z  =  R 

Hence  there  are  only  two  independent  variables,  since  R  is 
constant.  Alternately,  the  surface  may  be  described  by  two  inde¬ 
pendent  variables,  z  and  6  through  use  of  the  following  equations: 


3 


z 


same  z  as  before 


Here  z  and  0  would  be  called  the  "surface  coordinates."  Por 
other  examples  see  Figure  3.9  [ 1 3 . 


1.2  CLARIFICATION  OP  SELECTED  INPUT  DATA  CARDS 

This  Section  will  present  input  data  cards  which  were  judged 
to  require  some  additional  explanation  to  that  presented  in  the 
User's  Manual.  Not  all  cards  are  covered;  only  those  more  com¬ 
monly  used  or  specifically  requested  by  AFWAL/FIBR  are  discussed. 
If  the  explanation  in  the  User's  Manual  was  found  sufficient, 
then  the  card  will  not  be  discussed. 


Card  Variable  Comments 

( B-l )  IPOSTU  Only  used  for  plots,  data  written  to  FORTRAN 

unit  21 

/ 

IP0ST1  If  IP0ST1  =  0  saves  last  3  displacement 

vectors  on  unit  22 

If  IP0ST1  =  I  saves  last  I  displacement 
vectors  on  unit  22  (As  with  any  data  saved 
during  the  run,  the  user  Is  required  to 
furnish  the  proper  job  control  language  (JCL) 
in  order  to  save  the  specific  files) 

(B-2)  NIMPFS  Tape  or  unit  25  is  the  same  as  unit  22  saved 

from  a  previous  buckling  run.  It  is  recom¬ 
mended  that  Imperfection  amplitudes  be  on  the 
order  of  one  percent  of  the  panel  thickness. 
Note  that  only  local  w-displacements  for  each 
shell  unit  are  used  to  define  the  imperfection 
shape.  The  total  imperfection  is  printed  in 
the  STAG31  output. 

(B-3;  ..TAP  This  parameter  indicates  the  number  of  user 

parameters  required  by  the  user-written 
routines  if  used.  See  ( L— 1 )  card. 


(B-5)  WIMPPA 


There  are  NIMPPS  terms  on  this  record  which 
define  the  buckling  mode  amplitudes  used  for 
the  initial  imperfection.  Note  that  the 
buckling  modes  are  normalized  to  the  maximum 
value  in  the  eigenvector  and  usually  contain 
a  1.0  value. 

( C— 1 )  STLD(I)  The  base  load  is  defined  via  "Q"  cards.  Two 

different  base  load  systems  may  be  specified, 
load  system  A  and  load  system  B.  For  linear 
analysis,  the  total  load  on  the  structure  is 
given  by: 

Plin=  STLD(l)*  (Base  load  for  system  A)  + 
STLD(2)*  (Base  load  for  system  B) 

For  buckling  analysis,  load  system  B  is 
assumed  to  remain  constant.  The  total  load 
at  buckling  is  : 

Pbu  =  A*STLD(1 )* (Base  load  for  system  A)  + 

STLD (2 ) * (Base  load  for  system  B) 

Where  A  is  the  eigenvalue  found  from  the 
buckling  analysis.  For  nonlinear  analysis, 
first  a  linear  analysis  is  performed  using 
Plin  from  above.  Next  a  nonlinear  analysis 
is  performed  using  Pnin  from  above;  this  is 
called  step  1.  For  step  2,  a  nonlinear  analy¬ 
sis  is  performed  using  the  total  load  given  by: 

pnon=  (STLD ( 1 )+STEP ( 1 ) ) * ( Base  load  for  sys  A) 
( STLD ( 2 )+STEP ( 2 ) ) * (Base  load  for  sys  B) 

The  remaining  load  steps  would  then  normally 
use,  say,  (STLD(I)  +  STEP(I)  +  ...  STEP(I)) 
as  the  load  multiplier  unless  the  STEP(I)  is 
automatically  halved  or  doubled  by  STAGSC-1 
(see  Section  6  [1]  for  a  discussion  of  this). 
Here  1=1  for  load  system  A,  and  1=2  for 
load  system  B. 

If  the  postbuckling  response  is  being  calcu¬ 
lated,  the  starting  load  should  be  under  the 
linear  bifurcation  buckling  load  (e.g., 

STLD(l)  =  .5  *  Critical  Buckling  Load)  and 
the  load  step  size  such  that  2  or  3  steps  are 
required  to  reach  the  buckling  load  (e.g., 

STEP ( 1 )  =  .2  *  Critical  Buckling  Load). 
Additionally,  it  is  recommr  nded  that  STLD(I) 
and  STEP (I)  be  chosen  such  that  (STLD(I)  + 
STEP(I)  +...)  equal  about  1.0  at  the  estimated 
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( D— 1 )  ISTART 


NCUT 


critical  buckling  load.  This  is  also  the 
recommended  procedure  when  using  the  path 
length  as  an  independent  parameter  (nee  NSTRAT 
of  card  (D-l)).  Note  that  STEP(l)  and  STEP(2) 
have  to  be  greater  than  or  equal  to  zero. 
Negative  load  steps  are  not  allowed. 

For  a  new  case  or  for  using  the  POSTP  post¬ 
processor,  ISTART  must  be  set  to  zero. 

For  new  users,  a  value  of  four  is  recommended. 
As  the  user  gains  some  experience  and  eonfi- 
denece  in  nonlinear  analysis  this  value  may 
be  changed 


NEWT 


Same  comment  as  for  NCUT.  If  it  is  equal  to 
-20,  a  true  Newton-Raphson  method  is  used 
(update  the  stiffness  matrix  on  every 
iteration. ) 


NSTRAT  This  parameter  controls  the  solution  strategy 
used  in  the  analysis.  If  it  is  positive, 
then  a  modified  Newton-Raphson  method  is 
used.  If  it  is  equal  to  -1,  the  Riks  method 
will  be  used.  Traditionally,  the  postbuckling 
response  of  a  structure  has  been  obtained 
analytically  using  the  modified  Newton-Raphson 
solution  procedure  of  the  STAGSC-1  computer 
code  by  Incrementing  a  load  parameter  and 
solving  for  the  corresponding  solution  vector. 
Structural  instability  is  Indicated  by  a 
failure  of  the  modified  Newton-Raphson  solu¬ 
tion  procedure  to  converge,  even  with  a  very 
small  load  step,  or  by  a  change  in  sign  of 
the  determinant  of  the  tangent  stiffness 
matrix.  The  problems  associated  with  the  use 
of  a  load  parameter  as  an  independent  para¬ 
meter  are  related  to  solution  prediction 
beyond  a  limit  point  and  to  Ill-conditioning 
of  the  tangent  stiffness  matrix  in  the  neigh¬ 
borhood  of  a  limit  point.  The  Riks  method 
based  on  controlling  an  equilibrium-path-arc- 
length  parameter  Instead  of  the  traditional 
load  parameter  has  been  introduced  into  the 
STAGSC-1  computer  code.  The  Riks  method 
eliminates  one  singularity  in  the  tangent 
stiffness  matrix  at  critical  points  that 
causes  major  computational  difficulties. 

Using  the  Riks  method,  solutions  along  the 
unstable  equilibrium  path  of  the  postbuckling 
response  can  be  predicted  provided  only  one 
eigenvector  is  dominating  the  response  (i.e., 
there  is  no  modal  interaction  or  change  in 
buckle  pattern).  As  implemented  in  the 


STAGSC-1  computer  code,  the  Riks  method  may 
require  frequent  updating  and  refactorization 
of  the  tangent  stiffness  matrix  to  follow  the 
unstable  equilibrium  path  if  the  eigenvalues 
of  the  tangent  stiffness  matrix  are  closely 
spaced.  The  implementation  of  the  Riks 
method  is  such  that  the  load  step  size  will 
not  increase  after  it  has  been  cut  in  order 
to  trace  an  unstable  equilibrium  path.  The 
Riks  method  is  very  useful  but  requires  the 
analyst  to  carefully  watch  the  progress  of 
the  analysis. 

DELX  This  parameter  often  has  to  be  changed  on 

restarts  in  order  to  allow  a  nonlinear  solu¬ 
tion  to  proceed.  If  this  error  tolerance  is 
too  large,  and  the  stiffness  matrix  is  be¬ 
coming  nearly  singular  due  to  a  buckling 
phenomenon,  then  the  small  errors  accumulated 
tend  to  cause  numerical  problems  with  the 
stiffness  matrix  unless  DELX  is  lowered. 

WUND  The  use  of  this  factor  was  not  discussed  in 

the  User’s  Manual.  WUND  is  an  overrelaxation 
factor  that  the  program  changes  internally. 

It  is  recommended  that  the  program  be  allowed 
to  select  its  own  value  (i.e.,  do  not  input  a 
value  for  WUND). 

(D-2)  Note  that  if  IPRINT  is  set  to  a  value  of  1  or 

3,  the  buckling  modes  are  not  written  to 
TAPE22 . 

(G-2)  Partial  compatibility  contraints  are  defined 

via  a  base  freedom  and  a  slave  freedom  that 
is  equated  to  the  base  freedom.  These  con¬ 
straints  are  shown  in  the  computational 
degree-of-freedom  table.  Note  that  these 
constraints  apply  only  to  the  nodal  freedoms 
and  not  to  the  variation  between  nodes  (i.e., 
not  enforced  at  midside  nodal  d.o.f.). 

(G-4)  CC(I)  The  coefficients  used  for  the  constraints 

should  be  of  the  order  one  because  of  how  the 
terms  of  the  stiffness  matrix  are  scaled 
during  decomposition. 

(1-2)  Be  careful  in  defining  material  properties. 

The  nomenclature  used  in  STAGSC-1  is  not  the 
same  as  used  In  many  books  ~>n  composites.  A 
check  is  that  if  El  >  E2  then  V21  will  be  the 
larger  Poisson  ratio,  i.e.,  V21>V12. 


7 


(.1-1) 


SCY 


( K-l ) 
(M-2A ) 

(N) 

( Q— 3  )  LT 


The  "bar"  axis  used  in  defining  the  section 
properties  will  be  oriented  (using  the  "0" 
cards)  with  respect  to  the  "prime"  axis 
defined  at  each  row-column  intersection  or 
"node"  (See  Figure  3  [1]). 

For  an  example  of  calculating  these  factors 
see  Section  2.3.1. 

Note  that  a  maximum  value  of  NLAY  is  50. 

All  angles  are  in  degrees.  PRQP(5)  for 
ISHELL=5  (cylinder)  is  R  not  R2 . 

These  cards  enable  a  user  to  define  irregular 
grids,  specialized  grids,  or  uniform  grids. 

Irrelevant  for  initial  conditions  but  cannot 
be  zero. 


(R-l)  This  card  defines  most  of  the  output  informa 

tion  that  is  desired.  However,  additional 
stress  output  control  can  be  obtained  with 
cards  (K-2)  for  walls  and  some  of  the  stif¬ 
fener  or  ring  cards  (J-l). 


1.3  OUTPUT  AND  CLARIFICATION  OF  SELECTED  OUTPUT  INFORMATION 

This  section  discusses  interpretation  of  selected  output 
information  from  the  STAGSC-1  program.  Information  judged 
evident  or  obvious  has  been  omitted.  The  emphasis  here  is  to 
discuss  topics  that  the  user  of  other  finite  element  codes  may 
not  be  familiar  with  and  also  to  present  a  few  self-checks  for 
the  inexperienced  user. 

A  sample  input  job  control  language  (J CL)  to  run  the  VAX 
version  of  STAGSC-1  Is  presented  in  Figure  1.  The  use  of  various 
FORTRAN  units  is  Indicated.  The  preprocessing  program  which  must 
be  run  at  least  once  for  each  new  data  deck  is  STAGS1.  If  de¬ 
sired,  the  output  from  STAGS1  may  be  saved  from  FORTRAN  unit  2 
and  used  in  subsequent  restarts.  This  was  not  attempted  for  this 
report  since,  generally,  the  preprocessing  was  inexpensive. 

As  shown  in  Figure  1,  after  running  STAGS1,  STAGS2  is  run. 
3TAGS2  is  the  actual  finite  element  solution  procedure.  The 
output  from  each  program  is  now  discussed. 


«*  /  i 


$ 

Si  PROCEDURE  TU  EXECUTE  STAGS!  AND  STAGS2 

S 

S  SET  DEF  (.STAGS.) 

S  SET  VERIFY 
S 

$1  COMMENTS  -  JOBONE.INP  a  INPUT  DATA  FILE 
Si  JUflONE.RST  s  FILE  TO  SAVE  DISFS  FOR  RESTART 

SI  OUTPUT  WILL  BE  SENT  TO  SYSTEM  PRINTER 

S 

i  ASSIGN  J08UNE.INP  FOR005 
$  HUN  (STAGSCl)STAGSl.EXE 
S  ASSIGN  JOBUNt.RST  FQR02Q 
S  RUN  (STAGSCDSTAGS2.EXE 


Figure  1  Job  Control  Language  (JCL)  for  Executin'* 
STAGSC-1  on  the  VAX 


1.3.1  Preprocessing  (STAGS1)  Output 

This  output  Is  generally  not  dependent  upon  the  type  of 
analysis  to  be  performed  (linear,  buckling,  nonlinear,  etc.). 

Most  of  the  output  from  this  phase  of  solution  is  self- 
explanatory,  but  the  following  comments  are  included  for  the 
beginner : 

(1)  When  preparing  the  input  data  deck,  put  the  appropriate 
card  ldentif ication  ( ( B— 1 )  etc.)  as  a  comment  to  aid  in 
debugging  and  future  reference. 

(2)  Look  over  the  "INPUT  DATA  CARDS"  echo  for  any  obvious 
mistakes . 

(3)  Don't  worry  about  computer  memory  or  timing  information 
unless  problems  occur  or  timing  estimates  seem  very 
large. 

(4)  Check  the  mesh  SURFACE  COORDINATES  vs.  the  GLOBAL 
CARTESIAN  COORDINATES  to  see  if  they  make  sense. 

(5)  Check  ring  and  stringer  locations. 

(6)  Check  constitutive  matrices  to  see  if  they  appear  to  be 
the  right  order  of  magnitude. 

(7)  Check  loads  data. 

(8)  Note  the  COMPUTATIONAL  DOF  ASSIGNMENTS  for  future 
reference  (Section  1.3*2). 

Generally,  the  best  check  of  a  model  Is  to  run  a  simple  case 
where  expected  results  can  be  calculated  by  hand.  In  this  case 
the  output  from  STAGS2  (solution  phase)  might  be  checked  first. 
The  output  from  STAGS1  might  then  help  to  determine  possible  data 
errors.  However,  a  mistake  in  the  Input  will  sometimes  have 
little  effect  on  the  output  for  your  test  problem,  but  will  then 
cause  problems  when  the  real  problem  is  run.  Hence,  It  is  recom¬ 
mended  that  the  output  from  STAGS1  be  given  at  least  a  cursory 
check  before  the  model  is  considered  correct. 
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1.3.2  Solution  (STAGS2)  Output 

As  stated  above,  the  best  check  of  a  model  is  a  simple  test 
case  for  which  known  solutions  can  be  independently  calculated. 

For  complicated  structures  for  which  a  hand  calculation  may  be 
difficult  to  obtain,  there  is  at  least  one  test  case  that  can 
always  be  run. 

First,  a  set  of  boundary  conditions  is  picked  that  adequate¬ 
ly  prevents  rigid  body  motion  of  the  structure.  Next,  displace¬ 
ments  are  prescribed  for  those  boundary  points  (such  as  a  unit 
displacement).  The  output  should  show  that  all  displacements  of 
the  structure  are  equal  to  the  prescribed  boundary  displacements. 
This  may  be  generally  run  for  six  separate  load  cases,  each  load 
case  prescribing  a  unit  displacement  or  rotation  for  each  of  the 
coordinate  directions. 

With  reference  to  the  output  from  STAGS2,  the  following 
checks  and  comments  apply: 

(1)  Check  the  displacement  output  to  ensure  that  reasonable 
values  are  being  obtained  and  that  the  proper  boundary 
conditions  have  been  enforced.  The  displacements  for 
the  linear  solution  are  used  to  determine  initial 
stresses  for  a  buckling  analysis.  These  displacements 
are  also  used  as  initial  starting  vectors  in  a  nonlinear 
analysis . 

(2)  Equilibrium  forces  should  be  very  small  numbers  except 
for  points  where  boundary  conditions  or  applied  loads 
exist. 

(3)  Force  resultants  are  either  in  units  of  force/length  or 
moment/length. 

(4)  Stresses  in  composite  plates  depend  upon  the  ply  that 
is  being  examined.  Strains  often  produce  more  readily- 
interpreted  information. 


11 


(5)  The  MODE  DISPLACEMENTS  only  give  the  pattern  of  the 
buckled  structure  but  yield  no  information  as  to  the 
actual  displacements  that  occur  after  buckling;  a  non¬ 
linear  analysis  is  required  for  this.  The  modeshape  is 
usually  normalized  to  its  maximum  value.  Check  these 
displacement  patterns  to  see  if  reasonable  results  have 
been  obtained. 

(6)  In  nonlinear  analysis,  often  convergence  during  the 
first  few  steps  will  not  occur  if  the  first  load  step 
is  too  large.  If  this  is  suspected,  examine  the  step  0 
(linear)  displacements  and  the  step  1  (nonlinear)  dis¬ 
placements.  They  should  be  approximately  equal.  If 
they  are  not,  the  first  load  Increment  being  used  may 
be  too  large  and  cause  inaccurate  solutions  later  on. 
The  inaccurate  solutions  may  cause  nonconvergence  later 
in  the  analysis. 

(7)  Various  elements  have  additional  degrees  of  freedom 
associated  with  them.  These  additional  degrees  of 
freedom  are  discussed  briefly  in  Section  6  of  [1]  and 
in  slightly  more  detail  in  Section  4.3  of  [2].  The 
example  runs  in  this  case  used  the  411  element  which 
prints  out  an  additional  rotation  labeled  RW2.  This 
is  a  measure  of  the  Individual  elements  out-of-plane 
rotation  (see  Table  6.3.1  of  [1]). 

Diagnostic  messages  from  3TAGSC-1  are  generally  adequate 
for  debugging  a  model  although  often  one  has  to  hunt  through 
voluminous  output  to  find  them.  Section  7  of  [1]  should  be  read 
since  it  explains  most  of  the  non-obvious  messages.  One  comment 
which  applies  to  most  finite  element  codes,  including  STAGSC-1, 
will  now  be  discussed. 


Often  a  great  confusion  for  inexperienced  finite  element 
code  users  stems  from  messages  concerning  singular  or  nearly 
singular  stiffness  matrices.  One  output  from  STAGSC-1  is  the 
determinant  of  the  stiffness  matrix,  sometimes  accompanied  by 
comments  as  to  whether  the  matrix  is  ill-conditioned  or  not. 

The  determinant  of  the  stiffness  matrix  is  computed  as  follows. 
The  stiffness  matrix  is  decomposed  (factored)  into  the  form: 

K  =  UTD  U 


where 


U  is  an  upper  triangular  matrix  with  l's  on  the  diagonal. 

D  is  a  diagonal  matrix. 

Therefore,  det (K )=det (D ) ,  which  is  just  the  diagonal  terms  of  D 
multiplied  together.  Therefore  det(K)  depends  on  the  size  of  K 
and  will  be  very  different  if  a  100  x  100  matrix  or  a 
10,000  x  10,000  matrix  is  being  decomposed.  Therefore,  det(K) 
does  give  a  measure  of  the  singularity  of  K  but  not  really  a  very 
good  one.  Often  when  a  singularity  occurs,  STAGSC-1  also  prints 
out  the  smallest  diagonal  term  after  factorization.  This  is  the 
smallest  element  of  the  D  matrix  above.  By  noting  the  equation 
number  for  the  smallest  diagonal  term  and  cross-referencing  to 
the  COMPUTATIONAL  DOF  ASSIGNMENTS  (Section  1.3),  the  row-column 
location  of  a  potential  problem  area  may  be  found.  STAGSC-1  has 
an  option  to  remove  singularities  in  the  stiffness  matrix  which 
basically  constrains  that  degree  of  freedom  to  zero.  This  option 
is  appropriate  in  most  situations,  but  a  judgment  should  be  made 
to  ascertain  whether  the  other  displacements  and  EQUILIBRIUM 
FORCE  BALANCES  ((2)  above)  are  reasonable  before  accepting  the 
validity  of  the  results. 
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2.0  SAMPLE  PROBLEMS  INPUT,  OUTPUT  AND  INTERPRETATION 

In  this  section  various  sample  problems  are  presented.  The 
problems  are  presented  in  order  of  complexity.  The  problems 
include  flat  and  curved  plates,  stiffened  and  unstiffened,  linear 
static,  buckling  and  nonlinear  analysis. 

In  all  sample  problems  the  4ll  element  was  used.  Selecting 
the  best  element  often  depends  on  the  problem.  The  411  element 
was  chosen  here  based  partly  on  comments  in  Section  6.0  of  Ref¬ 
erence  [2]  and  partly  on  previous  experience.  Since  the  411 
element  has  been  found  to  perform  acceptably  for  most  problems, 
it  is  recommended  here  for  the  beginner.  Once  experience  has 
been  gained,  the  user  may  experiment  with  other  elements  in  order 
to  determine  their  relative  merits.  For  more  detailed  studies 
and  discussions  of  STAGSC-1  elements  the  reader  is  referred  to 
References  [2]  and  [6]. 

2.1  FLAT  COMPOSITE  PLATE 

This  section  describes  four  sample  runs  performed  on  the 
flat  square  plate  of  Figure  2.  The  dimensions,  layup  and  loading 
are  shown  in  Figure  2.  The  boundary  conditions  for  the  plate  are 
"classical"  simple  supports.  The  inplane  (X  and  Y)  translational 
boundary  conditions  are  shown  in  Figure  2.  The  input  and  output 
from  the  four  sample  runs  will  now  be  described. 

2.1.1  Input  for  Linear  Statics  (KNCOMPO) 

Figure  3  shows  the  input  for  a  linear  statics  run.  The  run 
is  given  the  label  KNCOMPO  for  future  reference  and  for  cross- 
referencing  to  the  computer  printout.  As  can  be  seen,  appropri¬ 
ate  comments  have  been  placed  In  the  Input  listing  of  Figure  3  so 
that  it  is  fairly  self-explanatory.  All  cards  in  the  input  have 
been  identified  as  to  their  type  using  the  STAGS  nomenclature 
(l'-w,  ( C— 1 )  etc.  For  a  full  explanation  of  any  particular  card 
type,  see  the  appropriate  data  card  (Section  3)  of  the  STAGS 
User's  Manual  [1], 
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Theoretical  Bucklinq  Load: 
N^r  =  15.66  Ib/in 
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LAYUP  =  [±45/90/0], 


STAGS : 


fTr  =  15.03  (KNC0MP1) 


=  19.6  x  10°  psi 
=  1.89  x  106  psi 
=  0.93  x  106  psi 
=  0.38  or  =  0.0366 
=  0.0056" 
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♦Buckling  of  Bars,  Plates  and  Shells,  Brush  &  Almroth,  p.  108, 
McGraw-Hill,  1975. 


Figure  2  Model  of  Flat  Plate  Used  in  KNCOMP  Series  of  Rt 
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add  in  (D-2)  and  (D-3)  cards  for  KNC0MP1  (Buckling  Analysis) 


The  first  card  must  always  be  a  comment  card.  The  second 
card  Is  called  the  Analysis  Type  Definition  Record  and  is  denoted 
as  card  ( B— 1 )  in  the  STAGS  User's  Manual.  The  "B"  type  cards 
define  overall  properties  of  the  model  such  as  how  many  different 
materials  are  being  used,  etc.  The  load  factor  of  1.0  on  the 
( G  —  1 )  card  is  the  factor  that  multiplies  the  load  system  that  is 
later  defined  on  the  "Q"  cards.  Material  properties  are  then 
defined  on  the  "I"  cards.  The  "K"  cards  define  the  number  of 
layers  and  the  layup  of  each  layer  for  each  wall.  They  are  re¬ 
peated  if  more  than  one  wall  ID  exists  (No.  of  shell  wall  is 
given  on  (B-3)).  The  "M"  cards  define  the  shell  type  (rectangu¬ 
lar  plate  in  this  example),  the  shell  dimensions,  and  the 
associated  wall  ID.  The  "N"  cards  define  what  element  type  is  to 
be  used  and  whether  reduced  or  standard  integration  is  to  be  used 
(reduced  integration  is  generally  recommended).  The  "P"  cards 
define  the  boundary  conditions;  the  "Q"  cards  define  the  loads; 
and  finally,  the  "R"  card  defines  output  options  and  control. 

The  changes  to  be  made  to  "KNC0MP0"  in  order  to  perform  the 
next  analysis,  which  is  a  buckling  analysis,  are  noted  in 
Figure  3.  Thus,  the  only  modification  required  in  the  present 
example  is  to  change  the  (B-l)  card  and  add  in  two  additional 
cards  (D-2)  and  ( D— 3 )  which  will  be  explained  under  the  linear 
buckling  input  description. 

2.1.2  Linear  Buckling  Analysis  (KNC0MP1) 

Figure  4  presents  the  input  for  a  linear  buckling  (or  linear 
bifurcation)  analysis.  The  linear  buckling  analysis  differs  from 
the  nonlinear  buckling  analysis  in  that  the  pre-stress  state  for 

the  linear  analysis  is  obtained  from  a  linear  stress  analysis. 

See  Reference  [1]  or  [3]  for  further  details  concerning  these 
differences.  As  can  be  seen,  most  of  the  input  cards  for  the 
buckling  analysis  are  the  same  as  the  linear  stress  analysis. 

The  main  difference  is  the  inclusion  of  the  (D-2)  and  (D-3)  cards 

and  the  first  entry  on  the  (B-l)  card.  An  abbreviated  explanation 

of  these  cards  is  included  in  Figure  .  Additional  information 
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may  be  found  in  Section  1.2  and  in  the  STAGS  User's  Manual.  The 
circled  items  in  Figure  4  indicate  modifications  to  be  made  for 
changing  to  a  nonlinear  analysis. 

2.1.3  Nonlinear  Analysis,  Triggering  Load  Included  (KNC0MP3) 

Figure  5  shows  the  input  for  a  nonlinear  static  analysis  of 

the  previously  described  plate.  Compared  to  Figure  4 ,  the  (B-l) 
card  is  changed  to  select  the  nonlinear  analysis  instead  of  the 
buckling  analysis.  The  (D-2)  and  ( D— 3 )  cards  have  been  removed 
since  they  only  apply  to  buckling  analysis,  and  the  D-l  card  was 
inserted.  The  (C-l)  card  has  been  modified  so  that  the  initial 
load  factor  is  1.0,  the  Increment  in  load  is  1.0  times  the  initial 
load,  and  the  maximum  load  factor  is  20.  The  (R-l)  card  was 
changed  although  this  was  not  necessary.  The  modification  of  the 
(R-l)  card  was  to  print  stress  resultants  every  other  step  as 
opposed  to  every  step.  Figure  5  also  includes  a  load  used  to 
trigger  the  initial  buckling  mode  of  the  flat  plate  being  analyzed. 
From  the  linear  buckling  analysis  (or  a  hand  analysis  in  this 
case),  the  initial  buckling  mode  may  be  determined.  Additional 
(Q-2)  and  (Q-3)  cards  were  then  added  to  apply  a  small  force 
(.001  lbs),  normal  to,  and  at  the  center  of  the  plate.  Also,  the 
(Q-l)  card  had  to  be  modified  to  allow  two  loading  systems,  and 
the  (C-l)  card  was  also  modified.  As  can  be  seen,  the  (C-l)  card 
applies  load  system  A  in  the  same  manner  as  before.  The  B  load 
system  however  is  "frozen"  at  its  initial  load  (max  load  factor 
on  (C-l)  for  load  system  B  =  1).  The  output  from  this  analysis, 
to  be  discussed  in  2.1.5,  shows  that  the  anticipated  result  was 
obtained . 

2.1.4  Nonlinear  Analysis,  Displacement-Controlled  (KNC0MP4) 

Figure  6  shows  the  Input  cards  which  now  allow  the  load  to 

be  applied  In  terms  of  end-shortening  or  end-displacement  of 
side  1  of  Figure  2.  As  can  be  seen  from  Figure  5,  the  only 
modification  required  was  that  of  the  first  (Q-^)  record.  The 
new  (Q-3)  now  specifies  an  end-displacement  of  5.0  x  10”-1  inches 
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as  the  initial  load.  The  reason  for  changing  to  a  displacement 
controlled  analysis  was  to  attempt  to  load  the  plate  through  the 
buckling  region.  In  general,  this  will  be  easier  using  displace¬ 
ment  control  since  the  stiffness  matrix  is  constrained  by  the 
specified  displacements;  thus  the  stiffness  matrix  should  not  be 
as  ill-conditioned  as  in  a  load-controlled  situation.  However, 
if  the  postbuckled  structure  softens  drastically  after  buckling, 
then  a  solution  through  postbuckling  may  be  quite  difficult  and 
require  many  restarts  to  accomplish.  The  present  example  is  such 
a  structure  and  will  be  discussed  further  in  Section  2.1.5. 

2.1.5  Discussion  of  Output  and  Results 

The  output  for  jobs  KNCOMPO  through  KNC0MP4  is  presented  in 
the  appendix.  Section  1.3  covered  general  interpretation  of  the 
output.  Section  4.0  contains  an  annotated  output  for  run  KNC0MP4, 
which  is  a  nonlinear  analysis.  Since  most  of  the  preprocessor 
(STAGS1)  output  and  much  of  the  solution  (STAGS2)  output  for  the 
linear  statics  (KNCOMPO)  and  linear  buckling  (KNC0MP1 )  are  the 
same,  the  reader  is  referred  to  Section  4.0  for  Interpretation  of 
this  output.  Additional  comments  on  the  buckling  output  (KNC0MP1) 
concerns  the  displacements.  STAGSC-1  prints  out  several  sets  of 
displacements  during  a  buckling  analysis.  The  first  set  labeled 
"DISPLACEMENTS  --  LINEAR  SOLUTION  — "  are  the  displacements 
obtained  using  a  linear  statics  solution  to  the  problem.  The 
next  set  labeled  "DISPLACEMENTS  MODE  i"  are  the  i * th  modes'  dis¬ 
placement  pattern,  and  are  not  actual  displacements.  They  are 
simply  the  displacement  pattern  and  are  usually  normalized  to  a 
maximum  value  of  1. 

Run  KNC0MP3  is  a  nonlinear  analysis  with  a  small  out-of- 
piane  load  (.001  lbs.)  used  to  trigger  the  buckling  mode.  Load 
control  was  used  to  increment  the  load.  The  plot  of  maximum  out- 
of-plane  displacement  versus  end-shortening  (Figure  7)  shows  that 
reasonable  results  are  being  obtained.  Figure  8  indicates  that 
the  load  in  the  skin  (Nx)  still  appears  to  be  linear  prior  to  the 
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theoretical  buckling  load  (15  lb/in).  If  the  analysis  were  con¬ 
tinued  through  buckling,  one  would  expect  Nx  to  begin  to  increase 
less  rapidly  with  increased  end-shortening.  Note  that  at  load 
step  9  (Nx  =  14.75)  the  solution  failed  to  converge.  This  is  due 
to  the  fact  that  at  the  buckling  load  (Nx  =  15),  the  stiffness 
matrix  is  singular  for  a  load-controlled  situation.  Figures  8 
and  9  were  plotted  using  a  separate  plot  routine  not  associated 
with  the  STAGSC-1  program. 

Run  KNC0MP4  was  an  attempt  to  use  displacement  control  on 
the  previous  problem  to  allow  solutions  to  be  obtained  through 
buckling.  When  a  structure  Is  loaded  using  displacement  control, 
the  stiffness  matrix  should  not  become  singular  for  this  problem, 
assuming  the  controlled  displacements  are  chosen  correctly.  This 
is  because  when  displacements  are  controlled,  the  stiffness  matrix 
has  the  controlled  degrees  of  freedom  eliminated;  thus,  if  the 
proper  displacements  are  controlled  (l.e.,  eliminated)  the  stiff¬ 
ness  matrix  will  be  rendered  nonsingular.  This  will  not  always 
happen;  for  example,  if  a  structure  has  several  buckling  modes 
that  are  not  properly  constrained,  then  the  stiffness  matrix  will 
become  singular  due  to  the  other  buckling  modes  forming.  In  this 
example,  the  results  are  actually  worse  than  those  using  load 
control.  This  is  therefore  another  example  of  the  often  quoted 
unpredictability  of  nonlinear  analysis.  The  displacement  con¬ 
trolled  approach  will  probably  yield  better  results  by  sharpening 
(decreasing)  the  convergence  criteria  (DELX)  and  forcing  the 
stiffness  matrix  to  be  updated  more  often.  This  approach  will  be 
discussed  in  Section  2.4  where  it  was  used  successfully.  Note 
that  the  results  obtained  for  KNC0MP3  and  KNC0MP4  are  not  "wrong;" 
the  solution,  however,  did  not  proceed  as  far  as  one  might  be 
interested  in  going.  The  Interested  reader  is  encouraged  to  run 
these  problems  and  attempt  to  "extend"  the  solution  obtained 
here. 


2.2  F! ,AT  COMPOSITE  PLATE  WITH  STIFFENERS 

This  section  describes  a  model  of  the  flat  plate  shown  in 
Figure  9.  The  plate  has  three  stiffeners  and  is  subjected  to  the 
shear  load  shown  in  Figure  9.  The  plate  dimensions,  layup,  and 
stiffener  properties  are  also  shown  in  Figure  9«  The  boundary 
conditions  will  be  described  in  Section  2.2.1.  This  sample 
demonstrates  one  option,  the  general  stiffener  section,  for 
modeling  stiffeners  and  also  demonstrates  the  use  of  different 
initial  and  incremental  boundary  conditions  on  a  structure.  Note 
that  the  use  of  either  the  general  stiffeners  or  of  the  subelernent 
stiffener  cards  precludes  certain  stiffness  distortions  (see 
Appendix  A  for  more  on  this). 

2.2.1  ohear  Buckling 

The  input  file  (CA2STIF)  for  the  present  problem  is  shown  in 
Figure  10.  Again,  the  Input  data  cards  have  been  defined  in 
terms  of  STAGS  nomenclature,  and  extensive  comments  are  present. 
Differences  between  this  Input  and  the  input  for  the  plate  of 
Section  2.1  will  be  described. 

The  Input  required  to  allow  stiffeners  (or  rings)  In  a  model 
affects  the  (B-2),  (B-3),  (P-2),  "J",  and  "0”  cards.  The  (B-2) 
and  (B-3)  cards  state  the  number  of  shell  units  with  discrete 
stiffeners  and  the  number  of  different  section  property  descrip¬ 
tions.  For  each  shell  unit,  the  (F-2)  card  states  how  many  rings 
or  stiffeners  are  present.  The  (J-l)  and  (J2-A)  cards  are  then 
repeated  for  each  section  property  description.  Finally,  the 
(0-2 A)  and  (0-2B)  cards  are  repeated  the  number  of  times  defined 
on  the  (F-2)  card;  i.e.,  for  each  stiffener. 

The  present  example  demonstrates  the  use  of  different 
boundary  conditions  for  the  "pre-critical"  stress  state  and  the 
"incremental"  solution.  STAGS  terminology  "pre-critical"  and 
"Incremental"  are  now  explained.  In  performing  a  linear  buckling 
analysis  two  matrices  are  needed,  the  small  displacement  stiff- 
matrix  K  and  the  so  called  geometric  stiffness  matrix  Kq. 
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K(j  is  formed  based  upon  a  set  of  stresses  and  is  a  linear  func¬ 
tion  of  the  stresses.  The  stresses  used  to  form  Kq  are  found 
from  a  linear  statics  solution  using  the  "pre-critical"  boundary 
conditions  specified  on  the  (P-1  and  2)  cards.  An  eigenvalue 
analysis  is  then  to  be  performed,  specifically  the  eigenvalue 
problem, 

( K  +  X  Kn  )  x  =  0 

U 

is  to  be  solved  for  a  specified  number  of  eigenvalues  X  and 
eigenvectors  (or  modeshapes)  x.  The  boundary  conditions  used  to 
form  the  matrices  in  the  above  equation  are  the  so  called 
"incremental"  boundary  conditions  and  are  specified  on  the  (P-3) 
card.  Appropriate  rows  and  columns  are  added  or  deleted  so  that 
Kq  and  K  are  the  same  dimension.  As  can  be  seen,  the  buckling 
analysis  is  "inconsistent"  in  that  the  stresses  are  calculated 
for  one  boundary  condition,  while  buckling  is  calculated  for 
another  boundary  condition.  However,  this  option  may  be  used  for 
modeling  more  complex  situations.  The  beginning  user  is  advised 
to  experiment  with  and  understand  this  option  before  utilizing 
it. 

2.2.2  Output  Discussion 

The  output  for  CA2STIP  is  contained  in  the  appendix.  The 
output  obtained  here  is  similar  to  the  output  discussed  previously 
in  section  2.1.5,  hence  discussion  will  be  brief.  One  point  to 
note  in  the  output  is  the  presence  of  a  negative  eigenvalue. 

Tnis  indicates  that  the  loads  applied  in  the  buckling  analysis 
must  actually  be  reversed  in  order  to  buckle  the  structure  in  the 
first  mode  (since  only  the  first  eigenvalue  is  negative  for  this 
o;t;ut j.  That  is,  if  the  loads  were  reversed  and  the  same 
nuew ling  analysis  were  repeated,  then  the  first  eigenvalue  would 
tr  positive  and  would  have  the  same  absolute  value  obtained  in 
to-'  appendix.  The  critical  buckling  load  in  the  present  case  is 
given  by: 


NCr  =  -3826.3511  x  1  lb/in  *  -3826.35  lb/in 


since  the  applied  load  was  1  lb/in  (see  (Q-3)  cards).  That  the 
applied  load  is  1  lb/in  may  be  verified  by  examining  the 
"EQUILIBRIUM  FORCES"  and  noting  that  the  accumulated  force  for 
row  1  in  the  Y  direction  is  20  lbs,  which  is  1  lb/in  for  a  20  in. 
wide  panel. 

2.3  CYLINDRICAL  SHELL  WITH  STIFFENERS 

Figure  11  presents  a  model  of  a  curved  panel  with  T  stiff¬ 
eners.  The  panel  dimensions  and  loading  are  shown  in  Figure  11. 
The  detail  of  the  stiffener  geometry  is  shown  in  Figure  12. 
Additionally,  the  auxiliary  "bar"  coordinate  system  which  is  used 
in  describing  the  stiffener  geometry  and  the  relation  to  the 
shell  "primed"  coordinate  system  is  shown.  This  example  illus¬ 
trates  the  use  of  a  cylindrical  surface  definition,  the  subelement 
option  for  defining  stiffener  and  ring  cross  sections,  and  the 
use  of  imposed  displacements  as  loading  conditions. 

2.3.1  Compression  Buckling  of  a  Cylindrical  Shell 

Figure  13  presents  the  appropriate  input  data  for  performing 
a  buckling  analysis  of  the  cylindrical  shell  of  Figure  11.  The 
shell  definition  is  affected  using  cards  (M-l),  (M-2A),  and  ( M — 5 ) - 
The  stiffeners  are  defined  as  in  Section  2.2,  except  (J-3B)  cards 
are  used  to  define  the  cross  section  instead  of  the  (J-2A)  cards 
used  in  Section  2.2.  The  (J-3B)  cards  are  repeated  for  each 
subelement  and  specify  the  locations  of  the  rectangular  subele¬ 
ments  in  terms  of  the  "bar"  axis  of  Figure  12.  Points  for  stress 
recovery  are  defined  on  the  second  (J-3B)  card.  The  "bar"  axis 
orientation  relative  to  the  "primed"  axis  (defined  at  each  node) 
is  specified  on  (0-2A)  cards  along  with  the  cross  section 
identification  numbers. 


node  on  shell 
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Figure  12  Detailed  Geometry  of  Stiffener 
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Figure  13  Input  for  CYL1  Problem 
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The  compression  load  is  applied  as  an  imposed  displacement 
of  .001  inches  on  side  1.  Additional  output  has  been  requested 
on  the  (R-l)  card  and  will  be  discussed  in  Section  2.3.2. 

2.3.2  Output  Descriptions  for  Compression  Buckling  of 
Cylindrical  Panel 

The  output  for  CYL1  is  contained  in  the  appendix.  Note,  in 
the  STAGS1  output,  the  SURFACE  COORDINATES  vs.  GLOBAL  CARTESIAN 
COORDINATES  (Page  6  of  output).  This  is  a  good  check  to  ensure 
that  the  shell  option  and  input  geometry  is  correct.  See  Section 

1.1.2  of  this  report  and  Section  2  of  [1]  for  further  discussion 
of  surface  coordinates.  In  order  to  validate  this  model,  an 
additional  run  with  a  small  modulus  input  for  the  stiffeners  was 
made.  The  results  of  this  analysis  (essentially  no  stiffeners) 
was  then  compared  to  the  exact  solution  for  a  curved  panel  [4J. 

2.4  ONE-BAY  MODEL  OF  STIFFENED  SKIN 

The  final  sample  problem  is  a  one-bay  model  of  a  stiffened 
skin  panel.  The  stiffened  panel  to  be  analyzed  is  shown  in 
Figure  14.  It  was  found  from  tests  of  the  panel  that  the  stiff¬ 
eners  were  so  stiff  that  buckling  of  the  skin  was  the  only 
concern.  The  buckle  pattern  was  established  and  it  was  found 
that  only  one-half  of  a  bay  needed  to  be  modeled  as  depicted  in 
Figure  14. 

This  sample  problem  illustrates  the  following  techniques  and 
capabilities  of  STAGSC-1  : 

•  Use  of  symmetry  in  modeling 

•  Use  of  linear  buckling  to  establish  preliminary 
buckling  load 

•  Use  of  imperfections  in  triggering  desired  buckling 
patterns  in  nonlinear  analysis 

•  The  restart  capabilities 


9*81 


iffened  Panels 


Since  this  report  is  aimed  at  the  beginner  in  finite  element 
analysis  some  comments  on  the  use  of  symmetry  are  in  order. 
Symmetry,  when  it  can  be  established,  enables  the  modeler  to 
model  only  a  portion  of  a  structure.  For  example,  in  the  linear 
static  analysis  of  a  simple  plate,  loaded  by  uniform  endloads 
along  the  length  of  the  plate,  only  a  quarter  of  the  plate  needs 
to  be  modeled.  Appropriate  symmetry  conditions  for  the  center 
lines  of  the  actual  plate  are  then  used  as  boundary  conditions 
for  the  quarter  model.  This  allows  not  only  a  much  more  economical 
solution  to  the  problem,  but  a  solution  which  is  actually  more 
accurate  than  if  the  entire  plate  were  modeled.  If  the  entire 
plate  were  modeled,  and  the  symmetry  conditions  were  not  obtained 
in  the  displacement  output,  the  differences  could  only  be  attrib¬ 
utable  to  numerical  round-off  or  the  inaccuracy  of  the  finite 
elements  used.  The  ability  to  detect  symmetry  conditions  will 
not  only  allow  the  engineer  to  perform  more  cost  effective 
analysis,  but  also  allows  the  detection  of  certain  modeling 
errors.  Additionally,  detection  of  symmetry  conditions  indicates 
good  engineering  Intuition  and  understanding  of  the  physical 
problem.  However,  in  buckling  analysis,  the  buckled  mode  must  be 
known  before  symmetry  can  be  used.  Erroneous  answers  will  result 
if  boundary  conditions  are  imposed  which  do  not  allow  the  correct 
buckle  pattern  for  the  entire  plate.  In  the  following  example, 
half  of  a  stiffened  skin  bay  was  modeled.  This  was  allowed  since 
the  buckled  pattern  was  known  prior  to  the  analysis. 

2.A.1  Linear  Compression  Buckling 

The  mathematical  model  is  shown  in  Figure  15.  Boundary 
conditions,  panel  geometry  and  ply  layups  are  also  presented. 

The  main  purpose  of  this  analysis  was  to  perform  a  nonlinear 
analysis  of  the  post-buckled  skin-stiffener  combination.  In 
order  to  get  a  feel  for  the  buckling  load,  a  linear  buckling 
analysis  was  first  performed.  The  input  (GDCPAN1)  is  presented 
in  Figure  16.  The  Input  is  well  documented  with  comments  and  is 
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Figure  15  Symmetric  One-Bay  Model  of  Stiffened  Skin 


very  similar  to  that  previously  discussed  in  Sections  2.1.2  and 
2.3.1.  The  output  for  the  buckling  analysis  is  labeled  GDCPAN1 
and  will  be  discussed  in  Section  2. A. 5. 

2.4.2  Nonlinear  Analysis  with  Initial  Imperfections 

As  discussed  in  [1],  it  is  often  necessary  to  use  initial 
imperfections  in  a  nonlinear  analysis  to  "trigger"  the  buckling 
mode.  This  was  found  necessary  for  the  present  panel.  The  first 
nonlinear  analysis  simply  loaded  up  to  several  times  the  buckling 
load  with  no  indication  of  nonlinearity.  Thus,  an  approximation 
to  the  buckling  pattern  shown  in  Figure  17,  from  GDCPAN1  was  used 
as  an  initial  imperfection  pattern.  Alternately,  the  buckling 
mode  calculated  in  Section  2.4.1  could  have  been  saved  on  TAPE22 
using  the  (B-2)  card  and  then  used  here  via  the  (B-2)  and  (B-5) 
card  combinations.  The  input  for  the  analysis  is  shown  in  Figure 
18  and  will  be  referred  to  as  GDCPAN7.  Besides  adding  in  the 
appropriate  cards  for  changing  from  a  buckling  analysis  to  a 
nonlinear  analysis  (discussed  in  Section  2.1.3),  the  following 
changes  are  made: 

•  Card  (B-l )  specifies  to  save  all  displacement  data 

•  Card  (C-l)  specifies  to  save  both  the  element 
(unfactored)  and  the  factored  stiffness  matrices 

•  Card  (M-6)  defines  the  initial  imperfections 

The  JCL  to  save  the  appropriate  files  for  this  analysis  is  shown 
in  Figure  19.  The  displacement  data  is  saved  on  GDCPAN7.S0D  and 
is  always  output  on  FORTRAN  unit  22.  Similarly,  the  element 
(unfactored)  stiffness  matrix  and  the  factored  stiffness  matrix 
are  output  on  FORTRAN  units  23  and  24,  respectively.  The  output 
for  this  analysis  will  be  discussed  in  Section  2.4.5. 
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i 

$1  PROCEDURE  TO  EXECUTE  STAGS1  AND  STAGS2 
$ 

S  SET  DEF  (.STAGS,) 

»  SET  VERIFY 

i  ASSIGN  GDCPAN7.1NP  F0R005 
S  ASSIGN  GDCPAN7.0U1  FOROOb 
5  RUN  (STAGSCDSTAGS1.EXE 
$  ASSIGN  GDCPAN7.0U2  F0R006 
i  ASSIGN  GDCPAN7 . SOD  FOR022 
i  ASSIGN  GDCPAN7 . KEL  FOR023 
S  ASSIGN  GDCPAN7.KFA  FQR024 
t  RUN  (STAbSCDSTAGS2.EXE 


ure  19 


JCL  for  Caving  Data  from  GDCPAN7 
for  Restart 


2. A. 3  Restart  of  Nonlinear  Solution  with  Stiffness  Matrices 
The  input  for  restarting  the  GDCPAN7  cun  is  presented  in 
'figure  20  and  is  labeled  GDCPAN8.  The  only  difference  between 
this  input  and  that  of  GDCPAN7  is  that  the  (C-l)  card  is  modified 
to  start  the  load  increment  at  one  of  the  previous  solutions 
(.928125  in  this  case),  and  the  stiffness  matrices  are  read  but 
not  saved.  It  is  noted  that  in  order  to  restart,  the  stiffness 
matrices  need  not  be  saved;  only  the  displacement  data  is  neces¬ 
sary.  Also,  the  (C-l)  card  specifies  that  a  load  step  .05  be 
used  compared  to  the  previous  value  of  .10.  The  (D-l)  card 
specifies  to  start  at  the  twelfth  step  of  the  previous  analysis 
(load  factor  =  .928125).  The  J CL  to  run  GDCPAN8  is  presented  in 
figure  21.  The  displacement  file  saved  as  GDCPAN7.S0D  is  equiva- 
lenoed  to  FORTRAN  unit  20  for  restart.  The  displacements  to  be 
saved  from  0DCPAN8  are  saved  on  GDCPAN8.S0D  and  are  equivalenced 
to  FORTRAN  unit  22.  The  stiffness  matrices  are  equivalenced  to 
the  same  units  that  were  used  to  save  them  (23  and  24).  Note 
that  STAGS  does  not  have  an  option  to  read  both  factored  and 
unfactored  stiffness  matrices  and  write  out  the  latest  factored 
and  unfactored  matrices  in  one  run. 

2.4.4  Restart  of  Nonlinear  Solution  without  Stiffness  Matrices 
The  input  presented  in  Figure  22  is  now  used  to  restart  run 
GDCPAN8.  The  (C-l)  record  indicates  that  the  initial  load  factor 
to  be  used  is  1.278125  and  that  stiffness  matrices  are  not  to  be 
read.  Thus,  only  the  displacements  will  be  used  In  the  restart 
procedure.  Also,  the  (C-l)  card  specifies  that  the  load  step  be 
.025  ani  that  loading  proceed  until  a  factor  of  1.75  is  obtained. 
Card  (D-l)  specifies  that  loading  will  start  with  the  eighteenth 
■•■ad  step  from  the  previous  solution  and  the  convergence  factor 
(iviLX)  has  been  reduced  to  1  x  10“^.  The  reduction  of  load  steps 
arid  allowable  convergence  error  are  typical  techniques  used  to 
improve  convergence  in  nonlinear  analysis.  Also,  the  updating  of 
the  stiffness  matrix  every  step  ((D-l)  card)  increases  the  rate 
of  convergence.  These  various  procedures  of  nonlinear  analysis 
are  discussed  in  Section  6  of  [1]  and  Chapter  6  of  [2]. 
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i 1  PROCEDURE  TO  EXECUTE  STAGS1  AND  STAGS2 
$ 

i  St  T  DEF  (.STAGS) 

$  SET  VERIFY 

i  ASSIGN  G0CPAN8.INP  FOROOS 
1  ASSIGN  GDCPAN8.0U1  FOR006 
l  RUN  (STAGSCDSTAGS1.EXE 
$  ASSIGN  GDCPAN8.0U2  FOR006 
i  ASSIGN  GDCPAN7 . SOD  FQR020 
i  ASSIGN  GDCPAN8 , SOD  FOR022 
i  ASSIGN  GOCPAN7.KEL  FOR02J 
S  ASSIGN  GDCPAN7.KFA  FOR024 
i  RUN  (STAGSCDSTAGS2.EXE 


Figure  21  JCL  to  Run  GDCPAN8 
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Figure  23  presents  the  JCL  to  run  GDCPAN9.  As  can  be  seen, 
the  displacements  saved  from  GDCPAN8  (GDCPAN8 .SOD)  are  equiva- 
lenced  to  FORTRAN  unit  20  while  the  current  displacements  are 
being  saved  on  FORTRAN  unit  22  in  case  another  restart  Is  required. 
The  output  for  this  run  will  be  discussed  in  Section  2.4.5. 

2.4.5  Output  Discussion  of  Half-Bay  Stiffened  Skin 

Runs  GDCPAN1,  8,  and  9  are  contained  in  the  appendix. 

GDCPAN1  is  a  standard  linear  buckling  analysis  which  predicts  a 
total  buckling  load  of:  P  =  AP'  =  4.584  x  158.15  -  725.  The 
value  of  P'  =  158.15  lbs  is  obtained  from  the  output  under 
"EQUILIBRIUM  FORCES  —  LINEAR  SOLUTION"  and  is  the  accumulated 
force  on  row  1  for  an  enforced  displacement  of  .001  in.  GDCPAN7 
is  a  displacement  controlled  run  which  uses  initial  imperfections 
obtained  from  GDCPAN1  to  trigger  the  buckling  pattern.  GDCPAN7 
was  run  until  a  load  step  of  .928  of  the ■ theoretical  buckling 
load.  At  this  point,  insufficient  computer  time  was  left  to 
continue  the  analysis.  The  stiffness  matrices  and  displacements 
vectors  were  saved  to  restart  the  problem. 

GDCPAN8  is  a  restart  of  GDCPAN7.  The  load  step  was  halved 
due  to  anticipated  problems  with  convergence.  The  run  continued 
loading  from  the  previous  .928  to  a  load  factor  of  1.28.  At  this 
point,  insufficient  time  was  again  the  reason  for  termination  of 
the  run. 

GDCPAN9  restarted  GDCPAN8  only  using  the  displacements  from 
the  previous  run.  The  load  was  successfully  applied  to  a  load 
factor  of  1.75  times  the  initial  buckling  load.  Normal  termina¬ 
tion  then  occurred.  Both  the  load  step  and  displacement  error 
(DKLX)  were  made  smaller  in  anticipation  of  convergence  problems. 

The  maximum  out-of-plane  displacement  vs.  end-shortening  is 
presented  in  Figure  24.  As  can  be  seen,  nonlinearities  occur  at 
an  end-shortening  of  about  .003  inches  which  is  less  than  the 
critical  end-shortening  (.005  in.)  predicted  using  the  buckling 
analysis.  Figure  25  shows  the  total  load  and  the  load  in  the 
stiffener  as  a  function  of  end-shortening.  As  can  be  seen,  a 
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si  PROCEDURE  TO  EXECUTE  STAGS1  AND  ST 

i 

S  SET  DEF  (.STAGS.) 

S  SET  VERIFY 

S  ASSIGN  GDCPAN9.INP  FOROOS 
$  ASSIGN  GDCPAN9.0UI  F0R006 
S  RUN  (STAGSCl)STAGSl.EXE 
S  ASSIGN  GDCPAN9 . 0U2  FOR006 
S  ASSIGN  GOCPANS.SGD  FOR020 
S  ASSIGN  GDCPAN9.SOO  FOR02J? 

S  RUN  (STAGSCDSTAGS2.EXE 


Figure  23  JOL  for  GDCPAN9 
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Maximum  Out-of-Plane  Displacement  for  GDCPAM 


slLght  decrease  in  linearity  and  load  occurs  near  the  end¬ 
shortening  of  .005  in.  which  is  as  anticipated.  The  stiffener 
load  remains  essentially  linear  due  to  its  high  stiffness  and 
constraints  imposed  along  its  length. 

For  comparison,  another  load  vs.  end-shortening  curve  for 
this  panel  is  presented  in  Figure  25.  This  curve  was  predicted 
using  the  PANCLP  [5]  program.  Comparing  these  two  curves  indi¬ 
cates  reasonable  agreement  between  STAGSC-1  and  PANCLP.  It  is 
noted  that  due  to  the  specialized  nature  of  the  PANCLP  program, 
its  associated  computer  cost  was  about  one-tenth  that  of  the 
STAGSC-1  program. 

3.0  ADVANCED  PROBLEMS 

Appendices  A  through  D  contain  four  examples  of  STAGSC-1 
advanced  modeling  capability.  These  examples  were  contributed 
by  Norman  F.  Knight  Jr.  of  NASA  Langley  Research  Center  whose 
helpful  comments  and  example  problems  are  here  acknowledged. 

The  sample  problems  are  considered  fairly  self-explanatory 
and  are  not  elaborated  upon  herein.  Appendix  A  is  a  study  of 
modeling  the  distortions  of  stiffeners  using  STAGSC-1.  Appendix 
B  is  an  example  of  how  to  model  a  multi-branched  structure. 
Appendix  C  is  an  example  of  modeling  a  complex  "pear"  shaped 
shell.  Appendix  D  presents  an  example  of  utilizing  the  "User 
Written  Subroutine"  capability  of  STAGSC-1.  These  examples 
should  aid  users  of  STAGSC-1  in  utilizing  these  more  advanced 
capabilities . 

4.0  ANNOTATED  OUTPUT  FOR  NONLINEAR  ANALYSIS  (KNC0MP4 ) 

Appendix  E  contains  portions  of  the  output  from  an  analysis 
previously  described  (KNC0MP4 ,  Section  2.1.4).  Since  it  is  often 
difficult  to  explain  something  without  seeing  it.  Appendix  E 
presents  the  output  with  abbreviated  explanations  printed  on  the 
actual  output.  Refer  back  to  Section  2.1.4  for  details  concerning 
the  model's  geometry,  loads,  etc. 
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APPENDIX  A 

REQUIRED  MODELING  DETAIL 
FOR  STIFFENED  PANELS 


REQUIRED  MODELING  DETAIL  IDENTIFIED  FOR  STIFFENED 
PANELS  WITH  RUCKLED  SKINS 


Norman  F.  Kni qht ,  Jr  . 
Structural  Mechanics  Branch,  SDD 
Extension  4585 
November  13,  1981 

(RTOP  505-33-33) 


The  desiqn  of  stiffened  panels  with  buckled  skins  is  of  considerable  interest 
to  the  aerospace  industry  because  of  the  significant  weight  saving  potential  ot 
buckled  skin  designs  when  compared  with  buckling  resistant  designs.  In  buckled 
skin  designs,  local  skin  bucklinq  is  allowed  which  results  in  additional  load 
being  transferred  to  the  stiffeners.  Because  postbuckling  deformation  shapes 
can  be  highly  complex,  the  analytical  modeling  detail  that  is  required  is  not 
known  in  advance  and  must  be  established  for  each  stiffened  panel  design.  One 
modeling  criterion  that  appears  to  have  merit  is  based  on  the  requirement  that 
the  model  used  to  predict  the  postbuckling  response  of  a  stiffened  panel  must 
first  predict  accurately  the  initial  buckling  response.  With  this  criterion  in 
mind,  a  study  was  made  to  identify  the  modeling  detail  required  by  a  nonlinear 
finite  element  code  called  STAGS  that  can  be  used  to  study  the  postbuckling 
response  of  stiffened  panels.  The  buckling  load  and  corresponding  mode  shape 
obtained  with  PASCO,  an  efficient  and  accurate  code  for  calculating  only  ini¬ 
tial  bucklinq  results,  were  used  as  the  standards  for  comparison  in  this  study. 
The  panel  selected  for  the  study  is  a  flat  graphite-epoxy  panel  with  a  16-ply 
skin  and  four  I-sti f feners .  The  buckled  cross  section  of  the  panel  as  deter¬ 
mined  by  PASCO  is  shown  in  the  upper  figure.  These  bucklinq  results  indicate 
that  the  stiffener  webs  deform  and  that  local  bending  occurs  near  the  skin- 
stiffener  interface. 

Several  STAGS  models  with  varying  levels  of  modelinq  sophistication  in  the 
skin-stiffener  region  are  shown  in  the  table  alona  with  the  buckled  cross  sec¬ 
tions.  In  the  analysis  of  stiffened  panels,  the  traditional  approach  has  been 
to  model  the  stiffeners  as  discrete  beams.  A  discrete  beam  model  accounts  for 
extensional,  bending  and  torsional  stiffnesses  and  any  eccentricities  of  the 
stiffener  by  lumping  these  properties  at  the  skin-stiffener  attachment  points. 
As  such,  the  discrete  beam  model,  STAGS-1,  neglects  not  only  the  cross  sec¬ 
tional  deformations  of  the  stiffener  but  also  any  local  bendinq  near  the  skin- 
stiffener  interface.  The  buckling  load  for  the  STAGS-1  model  is  10.6%  smaller 
than  the  PASCO  results  because  the  STAGS-1  model  has  a  larger  effective  stif¬ 
fener  spacing.  The  larqer  effective  stiffener  spacing  is  a  result  of  lumpinq 
the  properties  of  the  attachment  flanges  at  discrete  points.  In  order  to 
model  the  local  bendinq  near  the  skin-stiffener  interface,  a  plate  model, 
STAGS-?,  of  the  attachment  flanges  is  used  and  gives  a  bucklinq  load  which  is 
17. 9%  higher  than  the  PASCO  buckling  load.  The  bucklinq  load  is  hiqher  because 
the  discrete  beam  representation  of  the  stiffener  webs  used  in  the  STAGS-2 
model  does  not  allow  for  the  rolling  of  the  stiffeners.  The  next  refinement, 
STAGS-3,  treated  the  stiffener  webs,  in  addition  to  the  attachment  flanqes,  as 
plate  elements  which  then  allows  both  local  bending  and  cross  sectional  defor¬ 
mations  of  the  web  to  be  predicted  accurately. 

it  appears  that  the  level  of  modelinq  detail  required  for  accurate  analyses  of 
stiffened  composite  panels  which  are  designed  to  operate  with  buckled  skins 
is  greater  than  the  level  of  modeling  detail  that  previously  was  used  throuqhouf 
fhe  discipline  for  bucklinq  resistant  panels. 
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MODELING  OP  MULTI-BRANCHED  STRUCTURE 


5  BRANCH  STIFFENED  PANEL  MODEL 
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Sample  Problem  No.  1,  Pear  Shape  Cylinder 
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APPENDIX  D 


BRIEF  DESCRIPTION  OF 
USER-WRITTEN  SUBROUTINES 


In  this  Appendix  are  copies  of  routines  used  to  analyse  tie: 
postbuekiing  response  of  a  graphite-epoxy  curved  panel  witi 
a  central  circular  hole.  If  the  USRELT  routine  is  to.  or  us 
to  define  a  mesh  with  quadrilateral  elements,  then  the  foil 
ing  common  block  needs  to  be  included  if  the  default  values 
for  IN TEC  and  TPENL  are  not  used: 


COMMON/QUAFCM/  i  NTGP. ,  INTEC  ,  I  PENT  ,  TQRES  (  7  ) 


brief  Description  of  bscr-Wr  i  t  ten  Subi  emt  i  iu  s 

these  subroutines  are  used  to  model  the  square  pl_anform  oi  a  panel  witli  a 
single  central  circular  hole. 

i'he  finite  element  grid  is  generated  using  the  surface  coordinates  X  and  Y. 

The  grid  is  described  using  the  parameters  NKINCS  (the  number  of  rings  oi 
quadrilateral  elements  around  the  hole,  assuming  that  one  ring  will  be  for 
the  edge  supports  if  included)  and  NSPOKES  (the  number  of  radial  spokes  normal 
to  the  hole  boundary;  must  be  a  multiple  of  8). 


(X  ,  Y  ,  Z  )  .  .  Global  1:0- 
C  C  c  ,  .  , 

ordinates  of 

geometric  center 

of  panel . 

DHOLE  .  .  .  Diameter  of  hole 

RCURV  .  .  .  Radius  of  curvature 

of  panel  (set  equal 
to  zero  for  flat 
panel) 

A  ....  Length  of  panel 

EDGSUP  .  .  Location  of  edge 

supports  1 rom  edge 
of  panel 


NODE  NUMBERING 

The  node  numbering  scheme  is  counter-clockwise  beginning  at  the  edge  of  the' 
hole  as  shown: 


Two  additional  modeling  parameters  have  been  included:  MHKCR  and  KAT. 


MHRCK  .  .  .  Number  of  hole  radii  away  from  the  edge  of  the  hole  that  the 

rings  of  elements  will  be  circular. 

RAT  ....  Mesh  grading  factor.  Near  zero  gives  equal  spacing  of  rings 
increasing  the  value  of  RAT  to  a  maximum  of  one  gives  a  finer 
mesh  near  the  hole. 


BOUNDARY  CONDITIONS 

Internal  nodes  are  free.  Boundary  and  corner  nodes  read  in  as  data. 


I.OAPING 

Load  is  applied  using  a  point  force  at  the  top  (edge  1)  center  node  and 
applying  partial  compatability  constraints  across  edge  1. 

The  node  number  of  top  center  is  given  by 


NTC  =  NRINGS  *  NSPOKES  +  NSPOKES/2  +1 


Where  the  first  nodal  constraint  is  at  node  number 
PNC  =  NTC  -  NPC/2 

and  the  last  nodal  constraint  is  at  node  number 
LNC  =  NTC  +  NPC/2 

and  NPC  =  NSPOKES/4  (number  of  partial  compatability  constraints) 

REQUIRED  USER  DATA 

•  User  Parameter  -  Integers 
KQUAD ,  NRINGS ,  NSPOKES, 

IUVW  IRUVW  ]  repeat  for  each  edge,  edges  1,  2,  3,  4 
IUVW  IRUVW  ]  repeat  for  each  corner,  corners  1,  2,  3,  4 

IWALL ,  MHRCR 

•  User  Parameters  -  Floating  Point 

A,  DHOLE,  XC,  YC,  ZC,  RCURV,  RAT,  EDGSUP 

If  EDGSUP  =  0,  no  edge  support.  Along  the  lines  of  the  edge  support  W 
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LiiJtini,  u  T  Usur-Wri t ton  SubrouLiiit-s 


SUBKOUT  1  NE  IISKPT 
i; 

C  DEFINE  TUI'.  NODAI.  CDORI)  l  NATKS  AND  DOF 
C 

COMMON/  IIP!/  K(>UAI),  NR  I  NCS ,  NSPOKKS,  [BC(  I  f>)  ,  IWALL,  MHKCK 
DIMENSION  I  UVWB(  8)  ,  ll<UVWB(8)  ,DRAY(9A)  ,KRlNC(l)b) 

C 

C  KQIJAD  -  ELEMENT  TYPE :  IF  NEGATIVE,  USE  REDUCED  INTEGRATION 

C  NK1NGS  -  NUMBER  OF  INTERNAL  RINGS  OF  NOOKS  AROUND  THE  HOLE 
C  NSPOKKS  -  NUMBI.R  OF  RADIAL  SPOKKS  NORMAL  TO  UOLK  BOUNDARY 
C  (MUST  UK  A  MULTIPLE  OK  8) 

C  1UVWBCI-4)  -  U,V,W  BOUNDARY  CONDITIONS  ON  EDGES  OF  ELEMENT  UNIT 
C  IUVWB( 5-8)  -  U,V,W  BOUNDARY  CONDITIONS  ON  CORNERS  OF  ELEMENT  UN  LT 
C  IRUVWB( 1-4 )  -  kU.RV.RW  BOUNDARY  CONDITIONS  ON  EDGES  OF  ELEMENT  UNIT 
0  IRUVWB( 5-8)  -  RU.RV.RW  BOUNDARY  CONDITIONS  ON  CORNERS  OF  ELEMENT  UNIT 
C  I WALL  -  SELECT  SHELL  WALL  NUMBER  FROM  THOSE  GIVEN  BY  THE  K-CARDS 
C  MHRCR  -  MULTIPLE  OF  HOLE  RADII  FOR  CIRCULAR  RINGS  OF  ELEMENTS 
C 

COMMON/ P I E / DTR , RTD.PI 

COMMON/ UPF/ A , DHOLE , XC , YC , ZC , RCU  RV , RAT , EDCSU P 
C 

C  A  -  LENGTH  OR  ARC  LENGTH  OF  PANEL  ( SQUARE  PLANKORM) 

C  DHOLE  -  DIAMETER  OF  THE  HOLE 

C  XC.YC.ZC  -  GLOBAL  COORDINATES  OF  THE  CENTER  OF  THE  HOLE 
C  RCURV  -  RADIUS  OF  CURVATURE  OF  THE  PANEL 
C  ( . EQ.O  =>  FLAT  AND  .NE.O  =>  CURVED  ) 

C  RAT  -  MESH  GRADING  FACTOR 

C  (  NEAR  ZERO  GIVES  EQUAL  SPACING  OF  RINGS; 

C  AND  AS  RAT  =>  l,  FINER  MESH  NEAR  THE  HOLE) 

C  EDCSUP  -  DISTANCE  FROM  EDGE  OF  PANEL  TO  FIXTURE 
C  OK  BOUNDARY  CONDITIONS  OF  SUPPORT 

C 
C 

C  INITIALIZE 
C 

1 SY  S=  0 
IUVW=1 1 1 
IRUVW= 1 1 1 

IF(KQUAD.GE.420)  IRUVW--110 
POW=O.D 
IUS=0 
I  RS=0 
ICS=0 

RHOLE=0 . 5*  DHOLE 

CRAI) [=(  MHRCR+ 1  )*RMOLR 

DANG=2 .0*PI/ FLOAT( NS POKES) 

I  F(  EDCSUP.  LT. 0.(1)  F,nCSlJP=0.0 
DH=(  0  .  S*  A-KHOLE)  —  I’tM.SUP 
I  F( RAT. LE. 0.0)  RAT= 1 .OF-S 
IF(RAT.GT.l .0)  RAT=1.0 
RATE  I  =  RAT+ 1.0 

FACTOK=KAT/( P ATP  I  *  *  NR  I NCS- l .0) 


o  n  n  r>  n  n  r.  on  n  on 


P 14=0. 25* PI 

UNPACK  EDGE  AND  COKNKK  BOUNDARY  CONDITIONS  FROM  ARRAY  IRC 
JBC=  1 

DO  10  1=1,16,2 
IUVWB(  JBC)  =  l  BC(  I) 

10  JBC= JBC+ 1 

.!  BC=  l 

DO  11  1=2,16,2 
IRUVWB( JBC)= 1BC( I) 

11  JBC= JBC+ 1 

SET  UP  FOR  MESH  GRADING 

DR* DR* FACTOR 
DO  50  J=1 .NSPOKES 
ANG=( J-1)*DANG 
ANGR=ANG 

IF(ANG.GE.PI2)  ANGR=ANG-PI2 
IF( ANG.GE.PI)  ANGR=  ANG-P I 
IF(ANG.GE.(3.*PI2) )  ANGR=ANG-3.0*PI2 
IF( ANGR.LT.PI4)  BETA=( 1 .0/ C0S( ANGR)-1 .0) 

IF( ANGR.GE.PI4)  BETA=( 1 .0/ SIN( ANGR)-1 .0) 

DRAY( J)=( 1 .0+BETA) 

RRING( J)=RHOLE 
50  CONTINUE 

DEFINE  INTERNAL  NODES 

N  RI NGS 1=  NRI NGS- 1 

IF(EDGSUP.LE.O.O)  NRINGS1=NRINGS 
DO  100  1*1,  NRI NGS 
IM2= I- 2 

IF(  I.EQ.NRINGS. AND.EDGSUP.GT.0.0)  IUVW=110 

DO  100  J=l, NSPOKES 

IF( I.GT. 1 . AND.( RRING( l ) .GT. CRADI) ) 

.  RRING( J)=RRING(J)+DR*DRAY(J)*RATP1**IM2 
IF( I.GT. 1 . AND.( RRING( 1) .LE. CRADI) ) 

.  RRING( J)=RRING(J)+DR*RATP1**IM2 
RAY=RRING( J) 

IF( I.EQ.NRINGSI)  RAY=RRING( 1 ) *DRAY( J) 

IF( I.EQ.NRINGS. AND.EDGSUP.GT.0.0)  RAY=(0. 5*A-EDGSUP)*DRAY( J) 

GLOBAL  COORDINATES 

C  (X,Y)  ARE  SHELL  SURFACE  COORDINATES  AND  (XG.YG.ZG)  ARE  GLOBAL 
C  COORDINATES 
C 

ANG* ( J- 1 ) *  DANG 
X*  RAY*  CO  S( ANG ) 

Y*RAY*SIN(  ANG) 

XG*XC+X 


n  n 


IF(  RCIJRV.0T.O.O)  YOYC+Y 

1  K  RCUKV.GT.O.O)  /.(:-KCllKV*C()S(  AS  LN(  Y/KCUKV)  ) 
IK(RCURV.KQ.O.O)  YG-YO+Y 

LK(KOJRV.K<M).()) 

l  IJPT=(  I-  1  )*NSPOKKS+.l 

CALL  NOlUS(  1UPT,  IUS  ,  l  RS  ,  ICS  ,X(;  ,YG,  V.C. ,  1 UVW  ,  1  R1IVW  ,  LSYS.POW) 
10O  CUNT  l  Nil  K 
C 

C  DEFINE  Till;  IIOUNDARY  NOOKS 
C 

NS1=  3*  NSPOKKS/ 8+  1 
NS2=3*NSPOKKS/8H 
NS3R=NSPOKES/8+l 
NS3L=NSPOKKS 
N  S4=  7*  N  SPOK.ES/ 8+ 1 

<; 

C  SKT  IJP  EQUAL  LENGTH  SEGMENTS  ALONG  KOOKS 
C 

X=0.0 
Y=0 .0 

I) 0  200  1=  1  jNSPOKKS 

I I)  PT=  NRI NGS*  NSP(  >KE  S+ 1 
ANG=( I- 1)*DANG 
IF( I.LE.NS3R)  GO  TO  210 
IF( I.GT.NS3R.AND. I.LE.NS2)  GO  TO  220 
IF( I.GT.NS2.AN0.I.LK.NS1)  GO  TO  230 
IF( I .GT.NS 1 . AND. I .LE.NS4)  GO  TO  240 
IF( I.GT.NS4)  GO  TO  230 

RIGHT  SIDE  OF  EDGE  3 
C 

210  CONTINUE 
Il)CE=3 
I1JC07 
X=0.5*A 

Y=0. 5*A*TAN(  ANG) 

IF( I.EQ.NS3R)  GO  TO  310 
GO  TO  320 
C 

C  EDGE  2 
C 

220  CONTINUE 
I  KCE=  2 
I  BCC=6 

X=0.3*A*TAN( P 12- ANG) 

Y=0 . 5*  A 

I F(  I.KQ.NS2)  GO  TO  310 
GO  TO  320 
C 

G  EDGE  l 
G 

230  CONTINUE 


30 


I BCE* l 
[BC05 
X=-0 .  5*  A 

Y-0.5*A*TAN( Pl-ANC) 

IF( I.EQ.NS1)  GO  TO  310 
GO  TO  320 
C 

C  EDGE  4 
C 

240  CONTINUE 
I  BCE- 4 
IBCOB 

’  X=-0 . 5* A*TAN( 3.*PI2-ANG) 

Y=-0 . 5* A 

1F( I.EQ.NS4)  GO  TO  310 
GO  TO  320 
C 

C  LEFT  SIDE  OF  EDGE  3 
C 

250  CONTINUE 
IBCE=3 
X=0.5*A 

Y=-0 . 5*  A*  TAN( 2 .*  P I- ANG) 

C 

C  TRANSFORMATION  OF  COORDINATES  FROM  SHELL  SURFACE  COORDINATES 
C  TO  GLOBAL  COORDINATES 
C 

320  CONTINUE 
XG=XC+X 
YG=YC+Y 
ZG=ZC 

IF(RCURV.EQ.O.O)  GO  TO  321 
CC  YG=YC+RCURV* SIN( Y/RCURV) 

CC  ZG=RCURV*COS(  Y/RCURV) 

YG=YC+Y 

ZG=RCURV*COS(ASIN(  Y/RCURV)) 

C 

C  SET  UP  THE  EDGE  NODES 
C 

321  CALL  NODB( IUPT, IDS, IRS,ICS,XG,YG,ZG,IUVWB( IBCE) , IRUVWBC 1  BCE) , 
.ISYS.POW) 

GO  TO  200 
C 

C  TRANSFORMATION  OF  COORDINATES  FROM  SHELL  SURFACE  COORDINATES 
C  TO  GLOBAL  COORDINATES 
C 

310  CONTINUE 
XG^XC+X 
YG= YC+Y 
ZG=ZC 

IF( RCURV.EO.O.O)  GO  TO  311 


YC-YC+Y 


O-  kCUi:V*i  >>S(  AS  I  N(  V/  ItGUKV ) ) 

c 

c  si:  r  in’  tih:  corner  kih:k  nooks 
c 

i  1  1  CALI.  Noi)K(  IIIPT,  I  US,  LKS,  ICS,  XG ,  YG ,  ZC ,  IIIVWB(  I  BCG)  ,  I  KllVWItf  I  ICC)  , 

.  ISYS.POW) 

c 

200  CONTI  Nlli: 

c 

C  KNIT 
0 

RETURN 

END 

S  l  J  B  KOI  IT  INK  IISKK1.T 
C 

C  GENERATE  THK  MKSII  FOR  A  SQUARE  PLATE  WITH  A  CENTRAL  CIRCULAR  HOLE 

C 

COMMON/  UP  I  /  KQiJAl) ,  NK  (  NOS  ,  NSI’OKKS ,  lUVWB(H)  ,IKIIVWB(8)  ,  IWALL,  MI1RCR 

C 

C  KQ1IAI)  -  ELEMENT  TYPE:  IF  NECATIVE,  USE  REDUCED  INTEGRATION 
C  NRINGS  -  NUMBER  OF  INTERNAL  RINGS  OF  NODES  AROUND  THE  HOLE 
C  NS POKES  -  NUMBER  OF  RAD IAE  SPOKES  NORMAL  TO  HOLE  BOUNDARY 
G  (MUST  BE  A  MULTIPLE  OF  8) 

C  lUVWB(l-A)  -  U , V , W  BOUNDARY  CONDITIONS  ON  EDGES  OF  ELEMENT  UNIT 
C  IUVWB( 5-8)  -  U,V,W  BOUNDARY  CONDITIONS  ON  CORNERS  OF  ELEMENT  UN LT 
C  IRUVWB(  1-4)  -  RU.RV.RW  BOUNDARY  CONDITIONS  ON  EDGES  OF  ELEMENT  UNIT 
C  I RUVWB( 5-8)  -  RU , RV , RW  BOUNDARY  CONDITIONS  ON  CORNERS  OF  ELEMENT  UNIT 
C  IWALL  -  SELECT  SHELL  WALL  TYPE  FROM  THOSE  GIVEN  BY  THE  K-CARDS 
C  MHRCR  -  MULTIPLE  OF  ROLE  RADII  FOR  CIRCULAR  RINGS  OF  ELEMENTS 
C 

COMMON/ P I E / DTR , RTD , P I 

COMMON/ QUAFCM/ INTGR , I  NT KG , I  PEN L , IQKKS( 7) 

DIMENSION  ZKTA(96) 

G 

C  INITIALIZE 

i ; 

KCZ=0 .0 
IL1N=0 
IPLAS=0 
1NTEG=0 

I F( KQUAD.LT. 0)  INTKG= 1 
KQIJA1>=  I  ABS(  KQUAD) 

IPENL=0 

C 

C  ORIENTATION  ANGLE  IN  DEGREES 

C 

DANG=-2.0*PI*KTD/FLOAT(  NSPOKKS) 

Z.-/TA(  QO.O 
DO  ‘>0  I=2,NSPOKLS 
ZETA( I)=ZETA( I-l)+DANG 
‘)ll  CONTINUE 

C  FORM  THE  ELEMENT  GRID 

c 


8? 


00  200  L-  1  , NK  1  NOS 
L0DD=0 

IF(MOI)(  1,2)  .12).  I  )  IOUIM1 
DO  100  .1=1  .NSPOKES 
JODM) 

I  K(  MOI)(  J  , 2)  . K(>.  1 )  .!()l)l>=  1 
N  1  =  (  I-1)*NSP0KKS  +.( 

N2=ni+nsim>kks 

N3=N2+1 

N4=Nl+l 

1 F( J.EQ. NSPOKES)  N3=N!  +  1 
I  E(  J  .  EO. NSPOKES)  N4=N3— NSPOKES 
IF(  KQUAD.1.T.420)  00  TO  70 
IF(  1001). KQ.O.  AND. . 1001). EQ.  1)  GO  TO  7  5 
t F(  IODO.KQ.  1  .  AND. .JODI).  EQ.O)  00  TO  75 
70  CALL  QUA1KN1  , N2 , N3 , N4 ,KQUAD , IWALL,ZKTA( J)  ,KC 
.  ILIN, IPLAS) 

CO  TO  100 
75  CONTINUE 

ZET AR=  ZET  A( J ) - 90 . 0 

CALL  QUAD(  N2 ,N3 ,N4 ,N1  .KQUAD, IWALL.ZETAR ,ECZ  , 
.  ILIN, IPLAS) 

100  CONTINUE 
200  CONTINUE 
C 

C  EXIT 

C 

RETURN 

END 
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LisLin.;.s  oi  STACSC-1  input  data  and  S  IMM,  input  data . 


SQUARE  CY  I, I  NI)R  I  t'.AI.  PANEL  WITH  2-  I  NCH  CENTRAL  CIRCULAR  HOLE  '.A- I 
11  l  o  n  l  n  $8-1 

0  1  I)  0  ll)  0  4  All-2  ELEMENT  UNIT  ONLY 
2  0  4  29  $14-9 

0.00448  0.0896E-2  O.OH96E-3  0.0896E-4  $  14- $  I MPKKKECT I  <  )NS 


0. 

610192 

0.1 

2.0 

$C-1 

10  1 

l  $000 

20 

20  - 

■1  0 

.00001 

$o-i 

RIK'S  METHOD 

1 

421 

i) 

1 

1 

4  1 6 

0 

1 

$0-2 

U= CONSTANT  AT  X=0  CONSTRAINT 

1 

421 

0 

1 

l 

417 

0 

l 

$0-2 

1 

42  1 

0 

1 

l 

418 

(I 

l 

$0—2 

I 

421 

I. 

1 

1 

419 

0 

l 

$0-2 

1 

421 

0 

1 

1 

420 

0 

1 

$0—2 

1 

421 

0 

1 

1 

422 

0 

1 

$0-2 

1 

421 

0 

1 

1 

423 

0 

1 

$0-2 

1 

4  2  1 

0 

1 

1 

424 

0 

1 

$0-2 

1 

421 

0 

1 

l 

425 

0 

1 

$0—2 

1 

421 

0 

1 

l 

42b 

0 

1 

$0-2 

0 

0  0 

0 

0 

1  1 

$H- 

1 

ONLY 

IJS  KPT 

AND  USRELT 

1 

$1- 

l 

ISOTROPIC 

MATERIAL 

10 

.  Lb 

0.3 

$1- 

2 

2  $1-1  QUASI- ISOTROPIC  MATERIAL  (16  PLY  SKIN) 

19. ALA  0.0978  0.93K6  0.01  0.0  1.89E6  0.0  $1-2 

1  l  I  $K-1  SHELL  WALL  DLKI NITION 

1  0.08  0.0  $K-2 

2  1  L  A  $  K- I  SHELL  WALL 


•3 

L. 

0.005 

45.0 

$  K-2 

2 

0.005 

-4  5.0 

$K-2 

2 

0.005 

90.0 

$K-2 

r> 

0.005 

0.0 

$K-2 

2 

0.005 

45.0 

$  K-  2 

2 

0.005 

-4  5.0 

$K-2 

2 

0.005 

90.0 

$K-2 

2 

0.005 

0.0 

$K-2 

2 

0.005 

0.0 

$K-2 

2 

0.005 

90.0 

$K-2 

9 

L. 

0.005 

-45.0 

$K-2 

2 

0.005 

45.0 

$K-2 

2 

0.005 

0.0 

$K-2 

2 

0.005 

90.0 

$K-2 

') 

0.005 

-45.0 

$K-2 

0.005 

45.0 

$K-2 

4  1 

$K-1 

016=1)26 

0.08  $K-bA  WALL  THICKNESS 

O.A87AOAKA  0.2l8464Eb  0.  0.  0.  0.  ,  $K-5B 

o  .7 184A4EA  0.687606Kb  0.  0.  0.  0.  ,  S K—  S 14 

0.  0.  n . 2  44  A 7 1  LA  0.  0.  0.,  $K-5H 
o.  ...  ...  O.308650E3  0. I 4549K3  0.  ,$K-5B 
1 . .  o.  O.  0 . 14  5490L3  0 . 3607  28 K3  0.  ,  $K-$I4 

...  o.  0.  0.  0.  O.  1  $7 1  49L3  $K-  $  14 

4  1  1 A  K-  1  SHELL  WALL  (  +/-4S  ,90 ,0  ,0 ,90  ,-/  +  4  $)  SYMMETRIC 

2  0 .00  $6  4  $  .  0  $  K-  2 

2  0.0056  -4  $ .0  $K-2 

2  0 . 00  $ A  90.0  $K- 2 
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) 

o  .005t) 

0.0 

5K- 

’) 

0.0056 

0.0 

?K- 

2 

0.0056 

90.0 

$K-2 

2 

0.0056 

-45.0 

$K-2 

■) 

0.0056 

45.0 

$K2 

2 

0.0056 

45.0 

$K  2 

•) 

0.0056 

-45.0 

$K-2 

•> 

0.0056 

90.0 

?K-2 

•) 

0.0056 

0.0 

$E-2 

> 

0.0056 

0.0 

$K- 

2 

0.0056 

90.0 

$K-2 

2 

0.0056 

-45.0 

$K-2 

2 

0.0056 

45.0 

$K-2 

21  8  $  L- 1  NUPI+NUFF 

422,  SL-2A  UP!  VALUES 


10 

40, 

$  L-2A 

NRINGS 

AN1)  NSPOKES 

100 

000, 

$L-2A 

BC 

ON 

EDGE  1 

1  10 

100, 

$L-2A 

BC 

ON 

EDGE  2 

000 

000, 

$  L-  2  A 

BC 

ON 

EDGE  3 

1  10 

100, 

$L-2A 

BC 

ON 

EDGE  4 

100 

000, 

$  L-  2  A 

BC 

ON 

CORNER 

POINT 

1 

100 

000, 

$L-2A 

BC 

ON 

CORNER 

POINT 

2 

000 

000, 

$L-2A 

BC 

ON 

CORNER 

POINT 

3 

000 

000, 

$L-2A 

BC 

ON 

CORNER 

POINT 

4 

4,  SL-2A  SELECT  WALL  TYPE 

3  $L-2A  MULTIPLE  OF  HOLE  RADII  FROM  HOLE 

14.0  2.0  7.0  0.0  15.0  15.0  0.2  0.0  $L-2B  UPF  VALUES 

1  $U-1  LOADS 

1  1  0  $U-2 

2  5000.0  1  1  421  0  0  $U-3  POINT  LOAD  AT  NODE  421  X= 

1  0  0  0  0  0  $V-1  OUTPUT  FLAGS 

SQUARE  CYLINDRICAL  PANEL  WITH  2-INCH  CENTRAL  CIRCULAR  HOLE 
53  1  3.0  0  0  0  $PL-2 
0  $  P  L-  3 

0.25  0.0  0.0  0.0  0.25  $PL-4 

0  SPL-3 

0.25  -15.0  -45.0  -60.0  0.25  $PL-4 

2  10  0  1  0  0  1  0  0  0  $PL-3 


0. 

50 

90.0 

0 

.0 

0.0 

100 

.0 

10 

$PL-4 

2 

1 

0 

0 

1 

0 

0 

2 

0 

0 

-1 

$PL-3 

2 

1 

0 

0 

l 

0 

0 

3 

0 

0 

-1 

$PL-3 

2 

1 

0 

0 

1 

1 

0 

3 

0 

0 

-1 

$PL-3 

? 

l 

0 

0 

1 

5 

0 

3 

0 

0 

-1 

3  PL- 3 

2 

1 

0 

0 

1 

10 

0 

3 

0 

0 

-l 

$PL-3 

2 

1 

0 

0 

1 

15 

0 

3 

0 

0 

-1 

$PL-3 

2 

1 

0 

0 

1 

20 

0 

3 

0 

0 

-1 

$PL-3 

2 

1 

0 

0 

1 

25 

0 

3 

0 

0 

-1 

SPL-3 

2 

1 

0 

0 

1 

30 

0 

3 

0 

0 

-1 

3PL-3 

2 

1 

0 

0 

1 

35 

0 

3 

0 

0 

-1 

5  PL- 3 

2 

1 

0 

0 

1 

40 

0 

3 

0 

0 

-1 

3PL-3 

2 

1 

0 

0 

1 

4  5 

0 

3 

0 

0 

-1 

$PL-3 

2 

1 

0 

0 

1 

50 

0 

3 

0 

0 

-1 

$PL-3 

) 

1 

0 

0 

l 

55 

0 

3 

0 

0 

-1 

SPL-3 

0 

$  P  L-  l 
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AD-R153  142  SAMPLE  PROBLEMS  FOR  STAGSC-i<U>  ANAMET  LftBS  INC  SfiN 
CARLOS  CH  APPLIED  MECHANICS  DIV  P  J  UOYTOWITZ 
28  FEB  85  ANAMET-884.  1A  AFUAL-TR-84-2088 
UNCLASSIFIED  F33615-81-C-3201  F/G  9/2  NL 
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